The Ising ferromagnet in dimension five : link and spin overlaps 
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In the simple [hyper]cubic five dimension near neighbor interaction Ising ferromagnet, extensive 
simulation measurements are made of the link overlap and the spin overlap distributions. These 
"two replica" measurements are standard in the Spin Glass context but are not usually recorded in 
ferromagnet simulations. The moments and moment ratios of these distributions (the variance, the 
kurtosis and the skewness) show clear critical behaviors at the known ordering temperature of the 
ferromagnet. Analogous overlap effects can be expected quite generally in Ising ferromagnets in any 
dimension. The link overlap results in particular, with peaks at criticality in the kurtosis and the 
skewness, also have implications for Spin Glasses. 
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INTRODUCTION 

As is well known, the upper critical dimension (ucd) 
for Ising ferromagnets is four. Here we show results of 
simulations for an Ising ferromagnet with near neighbor 
interactions, on a simple [hyper]cubic lattice with peri- 
odic boundary conditions in dimension five, the next di- 
mension up. The Hamiltonian is as usual 
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with spins i and j near neighbors. Wc will take J 
and will quote inverse temperatures /3 = 1/T. 

All the principal properties of this system are well 
known. Recent consistent and precise estimates of 
the inverse ordering temperature j3 c are 0.11391 [lj, 
0.113925(12) 0, 0.1139139(5) from simulations, and 
0.113920(1) [41] from High Temperature Series Expansion 
(HTSE). We will use /3 C = 0.113915(1) as a compromise 
estimate. The susceptibility critical exponent and the ef- 
fective correlation length exponent take the exact mean 
field values 7 = 1 and v = 1/2, with a leading correction 
to scaling exponent 9 = 1/2 [5J. In the periodic boundary 
condition geometry, above the ucd the " effective length" 
is L cS = AL d l A where A is a non- universal constant 

We analyse simulation data for the "two replica" ob- 
servables link overlap and spin overlap, familiar in the 
Ising Spin Glass (ISG) context. It might seem curious to 
apply techniques developed for complex systems to the 
much simpler ferromagnet, particularly above the ucd, 
even though the observables can be defined in exactly 
the same way in a ferromagnet as in an ISG. Properties 
of the spin overlap at and beyond f3 c have already been 
studied in the 3d Ising ferromagnet Q. 

However, just because all the major parameters are 
well known it is convenient to use the 5d ferromagnet as 
a testbed for studying the critical behavior of various mo- 
ments or moment ratios (the variance, the kurtosis and 
the skewness) of both overlap distributions. The results 



should be generalizable mutatis mutandis to all Ising fer- 
romagnets. The final aim is to establish the ground rules 
for the properties related to the link overlap near critical- 
ity in order to apply a similar methodology to complex 
systems, in particular to ISGs. 

The link overlap 0] , like the more familiar spin overlap, 
is an important parameter in ISG numerical simulations. 
In both cases two replicas (copies) A and B of the same 
physical system, i.e. with identical sets of interactions 
between spins, are first generated and equilibrated; up- 
dating is then continued and the "overlaps" between the 
two replicas are recorded over long time intervals. The 
spin overlap at any instant t corresponds to the fraction 
q(t) of spins Si in A and B having the same orientation 
(Sf and Sf both up or both down), and the normal- 
ized overall distribution over time is written P(q). The 
link overlap corresponds to the fraction qe(t) of links (or 
bonds or edges) ij between spins which are either both 
satisfied or both dissatisfied in the two replicas; the nor- 
malized overall distribution over time is written Q(qe)- 
The explicit definitions arc 
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and 



*W = ^E^(i)S^(i)Sf (t)Sf (t) 
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where N is the number of spins in the system and Ni the 
number of links; spins i and j are linked, as denoted by 
ij. We will indicate means taken over time by (•••)• 

Equilibration and measurement runs were performed 
by standard heat bath updating on sites selected at ran- 
dom. The spin systems started with a random config- 
uration, i.e. at infinite temperature, and were grad- 
ually cooled and equilibrated until they reached their 
designated temperatures, where they received a longer 
equilibration. For example, for L = 10 the systems 
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at (3 = 0.114 saw roughly 10 6 sweeps before any mea- 
surements took place and the smaller systems at least 
10 7 sweeps. Several sweeps were made between mea- 
surements and with a flip rate of about 27% near f3 c 
which means at least four sweeps between measurements. 
For L = 10 about 10 7 measurements were recorded at 
each temperature near j3 c and considerably more for the 
smaller systems. 



LINK OVERLAP 

For any standard near neighbor Ising fcrromagnet with 
all interactions identical and with periodic boundary con- 
ditions, there is a simple rule for the mean link overlap 
in equilibrium (qe(f3,L)). If p s ((3,L) is the probability 
averaged over time that any given bond is satisfied, then 
by definition the mean energy per bond is 



\U(p,L)\ = l-2p a (p,L) 



(4) 



Because all bonds are equivalent 

{q t )(P,L) =p 2 s + (l-p s ) 2 -2 Ps (l- Ps ) = C/(/3,L) 2 . (5) 

This rule is exact at all temperatures (we have checked 
this numerically), so it would appear at first glance that 
link overlap measurements present no interest in a simple 
fcrromagnet as they contain no more information than 
the energy. However, the moments of the link overlap 
distribution reflect the structure of the temporary spin 
clusters which build up in the paramagnetic state be- 
fore /3 C , and the domain structure in the ferromagnetic 
state beyond f3 c . Thus if at some instant t a cluster of 
parallel spins exists in replica A and a similar cluster 
in the same part of space exists in replica B, then the 
instantaneous qe(t) will be significantly higher than the 
time average (qe(t)). The width of the overall distribu- 
tion Q(qi) increases rapidly on the approach to j3 c and 
we find phenomenologically that, as a consequence of the 
repeated occurrence of the cluster situation, around the 
critical temperature the distributions do not remain sim- 
ple Gaussians but develop excess kurtosis and skewness, 
even though to the naked eye these deviations from pure 
Gaussian distributions are not obvious; for instance no 
secondary peaks appear in the distributions. 

We exhibit in Figures 1 to 4 data at sizes L = 4,6,8 
and 10 for the Q- variance 
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the Q-kurtosis 
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FIG. 1: (Color online) The Q-variance, Eq. [S] as a function 
of size and inverse temperature for the 5d near neighbor fer- 
romagnet. In this and all the following figures the convention 
for indicating size is : L = 4, blue triangles; L = 6, red circles; 
L — 8, black squares; L — 10, pink inverted triangles. In this 
and following figures errors are smaller than the size of the 
points unless stated otherwise. The red vertical line indicates 
the inverse ordering temperature /3 C = 0.113915. 



and the Q-skewness 
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The three Q moments and moment ratios follow the 
standard definitions for the moments of a distribution. 
For the Q-variance we plot \og(Q vax ([3, L) — 1) so as to 
display the entire range of data. Fig. Q] shows the behav- 
ior of the Q-variance from high temperature through (3 C 
to j3 = 0.13, well into the ordered state. Fig.[2]shows the 
same data in the region near /3 C . 

It can be seen that the Q-variance has clear critical 
behavior. Just as for standard " phenomenological cou- 
plings" such as the Binder cumulant or the correlation 
length ratio £(/3, L)/L 5 / 4 (in the 5d case Q), it is size in- 
dependent at [3 C to within a finite size correction. As the 
Q-variance is a phenomenological coupling in this partic- 
ular system, it can be expected to have a similar form as a 
function of temperature in any fcrromagnet, with the ap- 
propriate finite size correction exponent. This makes the 
Q-variance a supplementary phenomenological coupling 
for ferromagnets in general. As qi is a near-neighbor mea- 
surement like the energy, the distribution Q(qe) tends to 
equilibrate faster on annealing than a parameter such as 
the correlation length £. 

The Q-kurtosis has a more unusual form, Fig. [3] At 
temperatures well above or well below the critical tem- 
perature it takes up the Gaussian value Qk(/3) = 3, but 
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FIG. 2: (Color online) The Q-variance as in Fig. [TJ in the 
region of the inverse ordering temperature /3 C . Sizes coded as 
in Fig. [T] 
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FIG. 3: (Color online) The Q-kurtosis, Eq. for the 5d near 
neighbor ferromagnet as a function of size and inverse tem- 
perature. Sizes coded as in Fig. [TJ 

near criticality there is an excess Q-kurtosis peak corre- 
sponding to a "fat tailed" form of the link overlap dis- 
tribution. With increasing L the width and strength of 
the peak decrease and the peak position /3 max ((9fc) ap- 
proaches (3 C . In the present 5d ferromagnet case the form 
of the temperature dependence evolves with L; with in- 
creasing L from a simple peak it tends to peak plus dip. 

The Q-skewness, Fig. 01 resembles the Q-kurtosis plot. 
The Q-skewncss starts at (a symmetric distribution) at 
(3 = and then develops a strong positive peak as a func- 



FIG. 4: (Color online) The Q-skewness, Eq. for the 5d 
near neighbor ferromagnet as a function of size and inverse 
temperature. Sizes coded as in Fig. [T] 

tion of j3 in the region of /3 C (so a distribution Q{qi) tilted 
towards high qg). Again the width and the strength of 
the peak decrease with increasing L. There is a weak in- 
dication of the beginning of a dip beyond the peak. The 
Q-kurtosis and Q-skewncss can be expected to show qual- 
itatively the same critical peak form in any ferromagnet. 

For the moment these observations are essentially phe- 
nomenological; it would be of interest to go beyond the 
argument given above in terms of correlated clusters of 
spins so as to obtain a full quantitative explanation for 
the details of the critical behavior of the Q(qe) distri- 
bution and its moments in finite L samples. Link over- 
lap moment peaks in ISGs resemble these ferromagnet 
results [13] implying that the peak structure is a very 
general qualitative form of the behavior of link overlap 
distributions at an Ising magnet critical point. 

The link overlap can be defined for vector spins fill ] by 

© = ^£[(Sf-S/)(Sf -Sf)] (9) 

ij 

which is invariant under global symmetry operations; the 
same link overlap critical properties as seen in Ising sys- 
tems may well exist in XY and Heisenberg magnets also. 



SPIN OVERLAP 

As for the link overlap, one can also define various mo- 
ments and moment ratios of the spin overlap distribution 
such as the P-variance 

P vai (P,L) = (q 2 ) (10) 
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FIG. 5: (Color online) The P-variance Eq.ITOlfor the 5d near 
neighbor ferromagnet as a function of size and inverse tem- 
perature. Sizes coded as in Fig. Q] 



before and at (3 C before taking on a two-peak structure 
in the ordered state [?|. This is in contrast to the mag- 
netization M-kurtosis (usually expressed as the Binder 
cumulant) in ferromagnets or the standard P-kurtosis in 
ISGs which both drop regularly from the Gaussian value 
3 towards the two peak value 1 with increasing order. 
However, it can be noted that the temperature variation 
of the kurtosis for the chiral order parameter in Hciscn- 
berg spin glasses has the same general form as the present 
ferromagnetic P-kurtosis, the distribution becoming fat 
tailed above the ordering temperature Q . 

As the distributions P{q) in equilibrium are by defini- 
tion symmetrical about q = 0, the P-skewness is always 
zero. However, for the one-sided distribution of the ab- 
solute value of \q\, other parameters can be defined, in 
particular the absolute P-kurtosis 

((M-(l<zl» 4 ) 
(M-(M))' 



and the absolute P-skcwncss 
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FIG. 6: (Color online) The P-kurtosis Eq.[TT]for the 5d near 
neighbor ferromagnet as a function of size and inverse tem- 
perature. Sizes coded as in Fig. Q] 
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The absolute P-kurtosis and the absolute P-skewness 
in the ferromagnet have rather complex phcnomcnologi- 
cal coupling temperature dependence patterns, with very 
weak finite size corrections at f3 c , Figs. [Jj and [8] If the 
weak finite size correction is a general property for these 
parameters, it could be usefully exploited so as to ob- 
tain high precision estimates of ordering temperatures in 
systems where these temperatures are not well known. 

The spin overlap properties do not transport from the 
ferromagnet into ISG systems in the same manner as the 
link overlap properties do because the P-variance takes 
on a different status : in an ISG q 2 becomes the order 
parameter. 



CONCLUSION 



and the P-kurtosis 

P fe (/3,i) = (g 4 )/< g 2 ) 2 (11) 

which is simply related to the Binder-like P-cumulant 
P b (p,L) = (3-P k (/3,L))/2. 

The P-variance in the ferromagnet has the phenomeno- 
logical coupling form, Fig. [5] at fi c the P-variance tends 
to an L independent value with a finite size correction. 
As the temperature goes to zero the P-variance will tend 
to L d . The P-kurtosis has a different phenomenological 
coupling form, a peak before a sharp drop to 1, corre- 
sponding to the P (q) distribution becoming "fat tailed" 



The standard near neighbor interaction Ising ferromag- 
net on a simple [hyper] cubic lattice in dimension five has 
been used as a test case in order to demonstrate the crit- 
ical form of the temperature variations of observables re- 
lated to the link overlap and the spin overlap, parameters 
developed in the ISG context and not generally recorded 
in ferromagnet simulations. The moments of the link and 
spin overlap distributions Q(qe) and P(q) show a rich va- 
riety of temperature variations, with specific critical be- 
haviors. The temperature dependence of the link over- 
lap kurtosis and the link overlap skewness show peaks at 
criticality which arc "evanescent" in the sense that they 
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FIG. 7: (Color online) The absolute P-kurtosis P\ q \, k , Eg. [HI 
for the 5d near neighbor ferromagnet as a function of size and 
inverse temperature. Sizes coded as in Fig. [T] 




will disappear in the the large size thermodynamic limit. 
The present results validate the assumption that these 
observables show true critical temperature dependencies. 
A temporary cluster phenomenon is proposed as deter- 
mining the evolution of the link overlap distributions. If 
correct, this mechanism is quite general, so we expect 
the critical peaks in the Q-kurtosis and the Q-skewness 
to be present in the entire class of Ising fcrromagnets, not 
only those in dimensions above the ucd, and plausibly in 
vector spin ferromagncts also. Beyond the class of fcr- 
romagnets, it can be noted that link overlap Q-kurtosis 
and Q-skewness critical peaks have also been observed 
in Ising Spin Glasses. 
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FIG. 8: (Color online) The absolute P-skewness P\ q \, s , Eq.[T3l 
for the 5d near neighbor ferromagnet as a function of size and 
inverse temperature. Sizes coded as in Fig. [T] 



